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The Monte Carlo program SUSYGEN, initially designed for computing distributions 
and generating events for supersymmetric particle production in e+e~ collisions, 
has now been upgraded to study supersymmetric processes at linear colliders by the 
inclusion of beamstrahlung, beam polarization, spin correlations and complex cou- 
plings including CP violating phases. It continues to offer, in the new context, the 
possibility to study the production and decay of supersymmetric particles within 
the most general minimal standard supersymmetric model (MSSM), the minimal 
supergravity model (MSUGRA) or the gauge mediated supersymmetry breaking 
model (GMSB), with or without assuming R-parity conservation. 



1 SUSYGEN2 content and upgrades 

Hopefully the linear collider would operate in an environment where a few 
sparticles would have been discovered and the high statistics will permit a pre- 
cise study of the soft supersymmetry breaking parametersEJ. One has therefore 
to move from a searches Monte Carlo to a "precision" physics Monte Carlo 
since one would need to measure accurately the above parameters. SUSYGEN2 
□ is a Monte Carlo event generator for the production and decay of supersym- 
metric particles and has been initially designed for e^e~ colliders. It has been 
extensively used by all four LEP experiments to simulate the expected signals. 
It includes pair production of charginos and neutralinos, scalar leptons and 
quarks. It offers also the possibility to study the production of a gravitino plus 
a neutralino within GMSB models and the production of single gauginos if one 
assumes R-Parity to be broken All important decay modes of supersymmetric 
particles relevant to LEP energies have been implemented, including cascades, 
radiative decays and R-Parity violating decays to standard model particles. 
The decay of particles is done through the exact matrix elements. The lightest 
supersymmetric particle (LSP) can either be the neutralino (xi), the sneutrino 
(j>) or the gravitino (G) in R-Parity conserving models, or there can be no LSP 
if R-Parity is violated. The initial state radiative corrections take account of 
Pt/pl effects in the Structure Function formalism. QED final state radiation 
is implemented using the PHOTOS □ library. An optimized hadronization in- 
terface to JETSElla is provided, which also takes into account the lifetime of 
sparticles. Finally, a widely used feature of SUSYGEN2 is the possibility to do 
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automatic scans on the parameter space through jiiscr friendly ntuples. Re- 
cently SUSYGEN2 has been upgraded to SUSYGENSB in order to adapt to the 
needs of the next generation of linear colliders, but also in order to extend its 
potential to supersymmetric particles searches at e^p colliders (e.g HERA) and 
hadronic colliders (e.g TEVATRON or LHC). The main new features relevant 
for linear colliders are the inclusion of beamstrahlung through an interface to 
CIRCEEEI, the full spin correlation in initial and final states, the inclusion of 
CP violating phases and the possibility to have an elaborate calculation of the 
MSUGRA spectrum through an interface to SUSPECT 



2 Beam polarization and spin correlations 

Since one expects high luminosities for the next generation of linear colliders 
(e.g. ~ 500/6~^ for the TESLA project), one can use beam polarization to 
reduce the standard model backgrounds and use the polarization dependence 
of the cross sections to study specific SUSY parameters. Moreover, as it has 
been stressed by several authors u, spin correlations play a major role in the 
kinematic distributions of final particles. To fulfill these two requirements, 
the helicity amplitude method a was used for the calculation of the different 
feynman amplitudes for production and decay, in order to obtain full spin 
correlation. Since such an amplitude involves products and contractions of 
fermionic currents, two basic functions, namely the B and Z functions were 
defined through: 

^Ai,A2(Pl'P2) = Ux,{pi,mi)PxUxJp2,m2), 

^XiM-MM^P^'P'^'P^'P*^ ^ [MAi(pi,mi)7''PAUA2(P2,TO2)] 

X [wA3(P3,"13)7^Pa'MAi(P4,TO4)] , 

where Pa stands for one of the two chiral projectors Pl or Pn and ux{p, m) de- 
notes the positive energy spinor solution of the Dirac equation for a particle of 
helicity A, four momentum p and mass to. The decomposition of the bispinors 
u\{p, to) in terms of the massless helicity eigcnstatcs LUx{k) yields analytic and 
easy to handle expressions for the B and Z functions. The amplitude is then 
factorized in terms of these basic building blocks; this fact permits compact 
and transparent coding and speed of calculation. The masses are not neglected 
in any stage of the calculation. For gaugino productions and decay, we use the 
widthless approximation. Let us illustrate it with, for instance, the calculation 
of the cross section associated to e+e~ X2X1 ~^ XiXi^'^'e". The total am- 
plitude associated to a given helicity configuration of the different particles is 
approximated by the product of the amplitude associated to the production of 
the two neutralinos (X2X1) with the amplitude corresponding to the decay of 
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the next to lightest ncutrahno (x2)- The propagator squared of X2 is approxi- 
mated by a factor given by Stt'^ / {rrij^oT ^0) . The total amplitude expression is 
given by: 

X^X?)xA^(x^^X?e-e+) 



The phase space integration is done through multichannel methodll3. The left 
side of figure illustrates the effect of spin-correlations in angular distributions 
of decay electrons for the neutralinos case^. Recent studiesEj have shown that 
T polarization effects yield valuable information for the MSSM parameters e.g 
for tan/3, the nature of the LSE, and the mixing angle Of. That is why, we 
interfaced SUSYGEN3 to TAUOLAllil which takes into account the r polarization 
accurately in the decays. The large effects one expects, for the momenta of 
the decay distributions are shown on the right side of figure |l|. 

M(i ) = 85 GeV/c^ M(x°) = 73 GeV/c', x", - B. taiip = 3 
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Figure 1: The figure on the left side shows the da /d cos 9 distribution for the final e~ asso- 
ciated to the process e+e~ — > X2X1 with and without taking into account spin correlations. 
The figure on the right shows the momentum distribution of e, /i and tt, decay products of r 
produced in the process e+e~ — > t^t^ 'r~^Xi''~~Xi- The distributions have been plotted 
for two hypotheses concerning the stau mixing angle Of (fj^ and Tfj) and compared to the 

no spin correlation case 
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3 Including phases in supersymmetry searches 

In the MSSM, there are new potential sources of CP non conservatioiEl. Com- 
plex CP violating phases can arise from several parameters present in the 
MSSM Lagrangian: the higgs mixing mass parameter fj,, the gauginos masses 
Mi , the trilinear couplings Ai . Experimental constraints on these CP violating 
phases come from the electric dipole moment of the electron and the neutron. 
Since all the couplings Ej and the different mass parameters fi, AIi, and the 
trilinear couplings A^, At and Ah have been assumed to be complex by default, 
the introduction of phases in the gaugino and sfermion sector for masses as 
well for cross sections has been straightforward in SUSYGEN3. Figure || shows 
the chargino pair production cross section variations in terms of 0^ , the phase 
associated to the fi parameter, for several values of the sneutrino mass rrij}^ 
and for a value of tan/3 equal to 1. One sees that there is a local minimum 
between the two extreme values = 0, 180 degrees), by default examined 
at LEP searches, but also that it is not so deep as to raise doubts on the 
exhaustiveness of "phaseless" searches. On the contrary at the linear collider, 
one cannot neglect this strong dependence of the cross section from phases. 
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Figure 2: e"'"e~ — > XiXi cross section evolution in terms of (f>^, the phase associated to 
the higgsino mixing mass parameter fi. We also assume different hypotheses concerning the 

sneutrino mass 



4 Future projects 

Native CH — h event generators are needed to optimize the couplings with linear 
collider simulation based on GEANT 4 and analysis tools. This will enable to 
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share code and reuse classes between different generators. The hehcity ampH- 
tude method and the two basic functions, namely the B and the Z functions, 
used to upgrade SUSYGEN seem to be optimally suited for object oriented imple- 
mentation. Moreover, Multichannel integration looks well adapted for object 
treatment. 
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